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With the goal of understanding the intricate behavior and dynamics of collections of neurons, 
we present superconducting circuits containing Josephson junctions that model biologically realistic 
neurons. These "Josephson junction neurons" reproduce many characteristic behaviors of biological 
neurons such as action potentials, refractory periods, and firing thresholds. They can be coupled 
together in ways that mimic electrical and chemical synapses. Using existing fabrication technolo- 
gies, large interconnected networks of Josephson junction neurons would operate fully in parallel. 
They would be orders of magnitude faster than both traditional computer simulations and biological 
neural networks. Josephson junction neurons provide a new tool for exploring long-term large-scale 
dynamics for networks of neurons. 
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How do large networks of neurons organize, communicate and collaborate to create the intrinsic behaviors and 
dynamics of the brain? Over the past century, individual neurons have been studied at the cellular, compartmental 
and molecular level. Synaptic models, while still somewhat rudimentary, accurately reflect many basic features of 
synapses and how they modify signals between neurons. Today, it is becoming feasible to explore networks of neurons 
behaving as units, how they synchronize, provide top-down or bottom-up feedback, and encode sensory information. 
This exploration is an important step toward understanding the brain, which will require multi-scale analysis with 
models of collective behavior at many different levels simultaneously. 

As part of this effort it is important that we understand how networks of neurons behave on the scale of thousands 
to tens of thousands of neurons, the size of a typical neocortical column. Large scale digital simulation projects 
such as the Blue Brain |2| and PetaVision [3[ projects have demonstrated that the limitations of inherently serial 
computer processors can be improved by effective parallel computing designs. But simulation time remains a significant 
hurdle to including biologically realistic features in large scale simulations. Analog simulations using VLSI circuitry 
to mimic neurons and synapses are improving in realism and speed, but they are still limited by complexity and 
power consumption. We propose a new direction for analog simulation of large scale networks of biologically realistic 
neurons. Using superconducting Josephson junctions to model neurons connected with real-time synaptic circuitry, 
we can explore neural network dynamics orders of magnitude faster than current digital or analog techniques allow. 
Using these circuits, we hope to learn about neural interactions such as synchronization, long term dynamics and 
bifurcations, feature identification and information processing. The long term goal is to understand group behavior 
of neurons sufficiently to use them as building blocks for studying larger scale neural networks and brain behavior. 

Our basic circuit unit (the JJ Neuron) involves two Josephson junctions connected in a loop as shown on the left side 
of Fig. [TJ The individual junctions behave phenomenologically like ion channels: one corresponds to a depolarizing 
current (like Na + ), and the other to a hyperpolarizing current (like K + ). Enhancements are possible, of course, 
and the inclusion of a third junction could allow for behaviors such as bursting that generally require at least three 
currents. The circuit displays many features of biologically realistic neurons such as the evocation of action potentials 
(firing) in response to input currents or pulses, input strength thresholds below which no action potential is evoked, 
and refractory periods after firing during which it is difficult to initiate another action potential 0, H| ■ 

The J J Neuron is a variation of well developed Rapid Single Flux Quantum (RSFQ) [y-Q circuitry and is thus 
straightforward to fabricate. RSFQ circuits using 20,000 junctions have been fabricated 0, so we estimate that a 
single chip could model as many as N — 10, 000 neurons, about the number of neurons in a cortical column. For 
larger brain regions, chips could be connected together. A simulated action potential takes about 50 picoseconds, 
and all neurons act in parallel. Fig [5] shows a comparison of the speed of a J J neuron with that of biological neurons 
and digital simulations using established models. The figure of merit displayed is the number of action potentials per 
neuron per second. Speed depends on the arithmetic complexity of the model employed, the number N of neurons 
simulated and their connectivity. We show speeds for sparse (no connections between neurons) and dense (all neurons 
connected to all others) networks as two extremes. While our estimates of speed are crude, they suggest that JJ 
Neurons are several orders of magnitude faster than either digital simulations or biological systems. This speed is due 
to three advantages: a) Josephson junctions have a shorter characteristic time than the conventional transistors used 
in computer processors; b) the JJ Neuron takes fewer cycles to simulate an action potential than digital simulations; 
and c) JJ Neurons compute in parallel using analog connections that do not affect computation time. 

This simple circuit simulates a single-compartment (space clamped) neuron and follows the spirit of mathematical 
neuron models such as Hodgkin and Huxley [Toj ] . Fitzhugh-Nagumo 

[llEl, Morris-Lecar pj], Hindmarsh-Rose 
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FIG. 1: Circuit diagram for the J J Neuron (left loop) connected to a model chemical synapse (right loop). In general, many 
synapses could connect to a single JJ Neuron. Orientation of junction phases is chosen for clockwise current in the left loop. 
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FIG. 2: Approximate number of action potentials(APs) per neuron per second for digital simulations, mammalian experiments 
and the JJ Neuron. FLOPs/AP is an estimate of the floating point operations required for one AP in each digital model and 
we assume a conservative 2 FLOPs per connection. We assume a CPU speed of 10 9 FLOPs/second for the digital simulations. 
The right three columns describe the speed of a network with N neurons. Sparse and dense connectivity represent extreme 
estimates of the computational time due to connections between neurons. Connections have no speed impact on experimental 
and JJ Neuron models as they are naturally parallel. 



and Izhikevich Each of these models describes the membrane potential (voltage) of the neuron as it varies in time. 
They also model, to varying degrees of detail, the ion channels which control the voltage via the flow of ions across 
the membrane. The ion channels open and close depending on the voltage thus providing the nonlinear feedback 
characteristic of neurons. Models differ in how many ion channels are described and their dynamic complexity. As 
with many of these mathematical models, our intent is not to model ion channel dynamics or chemical reactions in 
detail but rather to provide basic biological realism which allows us to explore the impact of network topology and 
connectivity strength. 

The JJ Neuron attempts to mimic important neuronal behaviors: action potentials, firing thresholds and refractory 
periods. The primary identifying feature of a neuron is the action potential (AP), i.e. a spike or pulse in voltage 
across the neuron membrane as the neuron "fires" . AP firing rates and inter-spike interval distributions are used to 
encode information in the brain [l5|. APs require at least two voltage-gated ionic currents, generally Na + and K + . 
The inward Na + current produces the rising phase of the pulse, while the outward K + current restores the membrane 



to its resting potential |16[. These currents have different time scales, resulting in the characteristic shape of the 



action potential. The dynamics of biological voltage-gated ionic currents are compl ex, and subsequent models have 
substantially expanded on those originally proposed by Hodgkin and Huxley [13, 03 ■ The J J Neuron, like most 
neuron models, simplifies the voltage dependence of ionic currents, modeling them by junctions whose dynamics are 
governed by familiar second order differential equations. It does include two ionic channels and supports APs. Other 
important features of a neuron model include the firing threshold and refractory period. The firing threshold is a 
level of external stimulus below which the response is negligible, while above it an AP is triggered. The refractory 
period is the period of time after triggering an AP during which it is difficult to evoke a second one. Biophysically, the 
refractory period is determined by the time it takes the ion channel proteins to "reset" to their initial configurations. 
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I. SINGLE NEURON MODEL 



Josephson junctions consist of two superconductors separated by a thin insulating barrier [19j. Electrons in each 
superconductor are described by a coherent wave-function with a definite phase. The phase difference from one side 
of the barrier to the other is the so-called Josephson phase 0, which controls all of the electrical properties of the 
junction, including the junction's voltage and current. Current can flow through the device without creating a voltage. 
This so-called super-current can be increased to the junction's critical current I above which a voltage develops given 
by V — (&o/2w)dcfi/dt, where t is time, $o = h/e is the flux quantum, h is Planck's constant and e is the electron 
charge. Normalizing current to this critical current, the current through a junction i depends on the phase <f> through 
the relation [2fjj : 

i = 4> + T<j) + sin(0). (1) 

The dot notation refers to differentiation with respect to normalized time r, where t 2 = i 2 $oC/27r/o, with C the 
capacitance of the junction. The normalized damping parameter is T 2 = $o/2ttIoR 2 C, where R is the resistance of 
the junction. The critical current, conductance (1/R) and capacitance of a junction are assumed to scale linearly with 
the cross-sectional area (A) of the junction. 

Equation [T] displays a useful analogy between the dynamics of a junction and that of a pendulum, allowing us to 
describe the action potential-like pulse generated by the circuit. Equation [1] is exactly the equation of motion for a 
damped and driven pendulum, where 4> represents the angle of deflection, T is the normalized drag constant, and i 
represents a driving torque. If i is constant in time, long term solutions include both static tilting (gravity balancing 
the torque) and whirling modes (where the torque overcomes gravitational forces). A third motion is also possible 
that combines these two. For appropriate parameters, and time dependent torque i, the pendulum will whirl over just 
once and settle back to a tilted state. In a Josephson junction, the phase whirling over just once creates a magnetic 
flux pulse, called a Single-Flux-Quantum pulse (SFQ pulse) 8]. This pulse has a similar shape and area each time it 
is stimulated. Pulses of this type form the action potentials in our neuron model. 

The JJ Neuron circuit shown in Fig. [TJ connects two Josephson junctions in a superconducting loop. This basic 
structure is a simplified superconducting digital component from RSFQ logic circuitry called a "DC-to-SFQ converter" 
Q. We call the two junctions the pulse junction and the control junction, denoted by subscripts p and c respectively. 
Two incoming currents (normalized to Iq p ) are called the input current ij n and the bias current i& which provides 
energy to the circuit. The signal to the synapse is the voltage across the pulse junction v p = 4> p , though the 
magnetic flux through the loop could be used for an inductive synapse circuit. The figure shows a model synaptic 
circuit connected above the pulse junction and thus driven by v p . The circuit parameters are the indicated branch 
inductances L s , L p , and L c which we scale by their sum L to tai to obtain A s , A p and A c . Using current conservation 
and fluxoid quantization [20j , we obtain two equations of motion: 

<P P + T0 P + sin(^p) = i p = -\(4> c + <P P ) + A s ii„ + (1 - A p )i b , (2) 
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+ T(j) c + sin(0 c ) =i c = -\((j> c + 4> p ) + Kkn ~ A p ib. (3) 



with coupling parameter A = $o/27rL tota ;/opi and geometric parameter 77 = A c /A p . To make the comparison with 
RSFQ circuitry more clear, DC-to-SFQ converters usually have A s « 1, A p = A c = and 77 < 1, whereas typical 
parameters for a JJ Neuron are A c = 0, A s = A p = 0.5, and 77 = 1. 

To explain how the circuit provides an action potential and then recharges to become ready to provide another, we 
describe the mechanism for generating a typical pulse. In equilibrium with ij n = 0, the bias current splits between 
pulse and control junctions and is large enough to put each junction just below its critical current. They are primed 
to achieve a whirling state. For ii n > 0, the input current acts to push the control junction away from whirling while 
pushing the pulse junction toward the whirling state. When the input current exceeds a threshold, it initiates the 
action potential: A voltage appears across the pulse junction creating magnetic flux in the loop which is analogous to 
a neuron's membrane potential rising. The flux in turn induces current in the loop pushing the control junction closer 
to its whirling state. As the flux builds, the control junction starts to whirl, draining the flux in the loop, ensuring 
that the pulse junction stops whirling and restoring the system so it can fire again. This process can repeat so long as 
the incoming current is held above threshold levels. The time lag between the pulse junction whirling, flux building 
and control junction responding to that flux is what creates the refractory period during which it is extremely difficult 
to initiate another pulse. The analogy goes beyond the membrane potential's relation to the flux $ = \(4> p + <p c ) 
normalized by L to tailop- The pulse voltage v p is analogous to a polarizing ionic current such as Na + while the control 
voltage v c is analogous to a hyperpolarizing ionic current such as K + . 
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The flux $ in the J J Neuron corresponds to the neuron membrane potential V m . Voltages across the pulse junction, 
v p , and control junction, v c , correspond to the Na + current I^ a and K + current Ik respectively. The input current 
tin corresponds to incoming post-synaptic current I syn . With these correspondences, the equation relating change in 
membrane potential to ionic currents can be written for the JJ Neuron 

~$ = v p + v c . (4) 

The analogous equation for a neuron is the Cable Equation for a space-clamped 2-channel neuron (l6j: 

CV m = Ik + In a + I syn ■ (5) 

The equations differ in that the synaptic current in ([5]), I sy m does not appear analogously as ii n in ([4]). Synaptic 
current does, however, appear implicitly since the phases <f> v and <fi c depend on through (O-©. The results of 
simulations below suggest that these terms have qualitatively similar effects. Like biological APs, the JJ Neuron AP 
is produced by the interaction of activating and restoring processes. 



II. SINGLE NEURON CHARACTERISTICS 



Single-neuron characteristics include the action potential, firing threshold and refractory period. JJ Neurons re- 
produce all three of these behaviors. Fig. [3] shows the voltage trace and ionic currents for APs in the J J Neuron and 
Hodgkin-Huxley models. 

Fig. [4] demonstrates the firing threshold. It shows the response to a brief current pulse of increasing strength along 
with an inset showing the threshold dependence of voltage peak on stimulus strength. Stimuli below threshold evoke 
a small (sub-threshold) response, while above the threshold full APs occur. 

Fig. [5] demonstrates the refractory period of the JJ Neuron. It shows the peak response times to two brief current 
stimulus pulses as the time between pulses varies. Current stimuli which are too closely spaced fail to produce a second 
AP. We emphasize that this refractory period is not explicitly hard-coded into the model (as for integrate-and-fire 
models [I]]]), but appears naturally as a result of differences in the time scales of the activating and restoring processes 
just as for biological neurons. The slight upturn in the curve for the second peak corresponds to interference between 
two closely spaced action potentials. In parts b-c are shown time profiles of APs responding to two brief current 
pulses. For large delay between pulses, the firing time of the second AP depends linearly on delay. When the delay 
is sufficiently short, no second AP occurs. 



III. CONNECTING NEURONS 



Synaptic models connect individual neuron models to form a network. Synapses control communication, signal 
transfers, and timing [22[. They can operate by means of direct electrical connections (electrical "gap junction" 
synapses) or chemical intermediates called neurotransmitters (chemical synapse). Electrical synapses are modeled 
by direct connection between J J Neurons. We focus on chemical synapses here as they are both more complex and 
more common. Chemical synapses can be excitatory or inhibitory, moving the post-synaptic neuron closer to or 
farther from threshold respectively. In chemical synapses, an AP causes the release of neurotransmitter molecules 
which diffuse across the synaptic cleft, bind to receptors on the post-synaptic membrane, and induce input currents 
in the post-synaptic neuron. The net effect on AP transmission across the synapse is a delay in transmission and the 
spreading of the pulse. We model the synaptic process using a resonant circuit attached to the pulse junction, shown 
on the right hand side of Fig. [TJ The output is taken across the capacitor and sent through a resistor to the input 
of a post-synaptic neuron. If the bias current applied to the JJ neuron is positive (negative) with respect to ground, 
then the synapse is excitatory (inhibitory). 

The equations describing our synaptic circuit and the coupling to the next neuron downstream come from Kirchhoff 's 
laws and involve normalized parameters describing the resonant frequency £Iq = t/ \JL S y n C syn , quality factor Q = 
QoRsynC syn / T i synaptic coupling A syn = L syn / L tota ^ and resistive coupling to the next neuron, r 12 , normalized by 
Ro p . We find that the voltage across the capacitor, v ou t (normalized by TIo p RO p ), and the current coupled to the 
next neuron, i\% (normalized by Iq p ) are given by: 
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FIG. 3: A. Time profile for action potentials in the JJ Neuron model. Parameter values for JJ Neuron calculations are A = 0.1, 
r = 1.5, A s = A;, = 0.5, 7] = 1. Input DC current is i; n = 0.21. B. Time profile for action potentials in the Hodgkin-Huxley 
model. The flux (black) in the JJ Neuron corresponds to the membrane potential (black) in the Hodgkin-Huxley model. 
Similarly, the voltages v p (red) and —v c (blue) in the JJ Neuron model correspond to the currents — Jjv a (red) and Ik (blue) 
in Hodgkin-Huxley model. Parameter values for all Hodgkin-Huxley calculations are C m = 1.01 uF/cm 2 ; G Na = 120 mS/cm 2 ; 
G K = 36 mS/cm 2 ; G L = 0.3 mS/cm 2 ; £ Na = 50 mV; E K = -77 mV; E L = -54.4mV; and T = 18.5°C unless stated otherwise. 
In this simulation an external stimulating current of 236nA is applied. The Hodgkin-Huxley neurons are modeled as isopotential 
cylinders with lengths and diameters of 500 um. 



"-syn \1 — ^syn ) : , ' 12 . . , \ ,„n 
^ 112 + — «12 = Vout - J^syn [V c2 + V p2 ) . (7) 

Here v p2 and v c2 represent the voltage across the pulse and control junction of the next neuron, respectively. 

In addition, the synapse circuit's back-action on the JJ Neuron circuit creates two additional terms on the right 
side of —i\ 2 — Xv ou t/ (A^^g). For the simple demonstrations of coupling behavior shown here, this back-action 
was not strong enough to warrant buffer junctions, as is commonly needed with RSFQ circuits [23| . However, for 
more complex circuits buffer junctions will most likely be needed. 

We demonstrate synaptic behavior for inhibitory and excitatory coupling in a two-neuron setting. The first neuron 
is connected to the second via a synaptic connection as shown in Fig. [TJ The results are shown in Fig. [B][7J 

In Fig. [6l we show excitatory coupling from neuron 1 to neuron 2. Neuron 1 receives an external stimulus, and its 
output drives neuron 2, which responds by firing at the same rate but out of phase with neuron 1. 

In Fig. we show inhibitory coupling from neuron 1 to neuron 2. Neuron 2 is configured to fire repetitively by 
increasing its bias current. When neuron 1 is not firing, this external stimulus causes neuron 2 to fire repeatedly. 
However, when neuron 1 is stimulated, it inhibits the firing of neuron 2 even though neuron 2 continues to receive the 
external stimulus. Both return to their original state after the stimulus to neuron 1 is removed. Qualitatively similar 
results (Fig. [5][5]) are obtained with two coupled Hodgkin-Huxley neurons. The back-action of neuron 2 on neuron 1 
causes small subthreshold oscillations and a slight decrease in the firing rate of the first neuron (~ 10%). These minor 
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FIG. 4: The threshold response of flux time-traces colored by the strength of impulse input signals. Stronger inputs lead to 
action potentials. The inset graph shows how peak response depends on the strength of the input pulse. Parameter values as 
for Fig. [3K except V = 1.0. Input current is a single square pulse of width 5 and varied height initiated at t = 30. 



effects are not present in the Hodgkin- Huxley model because the alpha function synapse does not allow back-action. 
Adding buffer junctions between the JJ neurons would be needed to obtain precise agreement. 

IV. DISCUSSION 

We have shown that the J J Neuron is a biologically realistic model for single neuron dynamics, and that individual 
J J Neurons can be coupled together in ways that are analogous to inhibitory or excitatory chemical synaptic coupling. 
We envision that large scale long term dynamics of networks of neurons can be explored using J J Neurons, contributing 
to our understanding of behaviors such as synchronization, pattern recognition, and memory formation. JJ Neurons 
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FIG. 5: a) Response to twin pulse inputs with delay shows the absolute refractory period. On the left, the time of the first 
(blue) and second (red) response peak is graphed against the delay between input pulses. The vertical line represents the delay 
below which no second pulse is created. On the right, time profiles are shown for input currents (green) and output voltages 
(blue). In b) two pulses are generated while in c), the delay is too small and no second pulse arises. Parameter values as for 
Fig. EJ\ except F — 1.0. Input pulses are twin square pulses of height 0.54, width 5 and varied delay between pulses. The first 
pulse initiates at t = 30. 
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FIG. 6: Excitatory synaptic coupling. One JJ neuron (blue) drives another (green). The system is initially quiet. At time 
50, a constant input current (red) causes the first neuron to fire repetitively. Excitatory coupling induces repetitive firing of 
the second neuron. Parameter values as for Fig. [3j\ except if, = 2.0, F = 2.0. Synaptic parameters are Q = 1 .0, Q = 0.05, 
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FIG. 7: Inhibitory synaptic coupling. The second neuron (green) is configured to fire repetitively. When the first neuron 
(blue) is activated by an external input current (red), it inhibits the second. Both return to their original state after the 
external stimulus to the first neuron is removed. Parameter values as for Fig. |6]except in = —1.9, ib2 = 1.76, A S 2 = 0.65, and 
A P 2 = 0.35. This pus the second neuron into a repetitive firing mode. Synaptic parameters are also the same except ri2 = 0.6. 
Input current to first neuron is a single square pulse of height 0.3, width 360, at t = 240. 



should be able to alleviate computational bottlenecks currently slowing needed simulations. We now discuss the 
advantages and limitations of this approach. 

The advantages of large scale JJ Neuron simulations are speed, biological realism, simplicity of circuit design and 
low power consumption. The major advantage is speed. Based on similar circuits constructed for RSFQ circuits @, 
network simulations of 20,000 densely coupled neurons are reasonable and could simulate one trillion APs for each of 
these neurons in a few minutes. This speed is unachievable with digital simulations using modern computers. Much 
of this advantage is due to our analog rather than digital approach. 

Analog simulations using electrical circuits to model neurons (24T - [30l | provide a framework more in line with the 
parallel nature of biological neural networks. Silicon-based Very-Large-Scale-Integrated (VLSI) circuits have been 
designed which simulate simple integrate and fire neurons (3lL |32| , bursting neurons (33j , and plastic synapses with 
timing and homeostasis [3lll34|. The focus of VLSI research has been on implementation of a learning processor rather 
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FIG. 8: Behavior of two Hodgkin-Huxley model neurons coupled by a simple excitatory synapse model. The synapse uses 
a positive "alpha" function [la ] with maximum current 400 nA and time scale 0.1 ms which is triggered at -30 mV on the 
downward slope of the presynaptic neuron's action potential. The presynaptic neuron (blue) is made to fire repetitively by 
means of a constant current input (red curve) initiated at time 5 ms. Excitatory coupling then causes repetitive firing of the 
postsynaptic neuron (green). 
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FIG. 9: Behavior of two Hodgkin-Huxley neurons coupled with an inhibitory synapse model using a negative alpha function 
with maximum current 215 nA and time scale 1 ms. The second neuron (green) has E L set to -15 mV to make it fire repetitively 
in the absence of external inputs, while the first has E L set at the standard -54.4 mV value to make it normally quiescent. 
When the first neuron (blue) fires as the result of an external input current (red), it inhibits the second. Both return to their 
original state after the external stimulus to the first neuron is removed. 



than simulation of long-term neuron dynamics so speed comparisons are difficult. Speeds of VLSI analog circuits are 
typically chosen to emulate biological neuron timescales (order milliseconds). Presumably these circuit speeds could 
be increased up to a few GHz. In contrast, RSFQ processors similar to the J J Neuron have been clocked on the order 
of 100 GHz Q. In addition, the power dissipation for a J J Neuron network should be significantly less than a VLSI 
circuit of the same size, allowing larger networks. Finally, the non-linear behavior of Josephson junctions allows a JJ 
Neuron to contain only two Josephson junctions, much smaller than the 22 transistors for a VLSI Integrate and Fire 
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neuron (3lj . 

The major limitation of JJ Neuron models compared to digital simulations is measurement. While every neuron in 
a computer simulation can be separately monitored, this is not technically feasible in large networks of JJ neurons. 
So for settings which require detailed monitoring of each cell in the network, the JJ Neuron is not appropriate. But 
many settings do not require such detailed observation, instead measuring averaged output or the output of a few key 
neurons. The measurable data from JJ Neuron experiments aligns more closely with measurements of living tissue. 
Placing magnetic field detectors on the chip allows measurement of averaged activity levels, and subsequent comparison 
of simulation data with electroencephalogram (EEG) measurements. Other settings which involve output from only 
a few key neurons are ripe for JJ Neuron models. Quantities such as voltage, current or flux at important points 
in the network can be measured directly, limited by the number of wires passing through the cooling apparatus. In 
addition, a variety of collective behaviors can be measured using already developed RSFQ logic circuits such as clocks, 
counters, splitters, followers, flip-flops and transmission lines [7|. So, while J J Neuron simulations are not well suited 
to some questions, they are for others, especially those involving averaged or collective behaviors or producing output 
from a small fraction of neurons. Analog simulations will not replace digital simulations. Instead, they complement 
each other, with analog simulations providing collective measurement on long time scales and large networks, while 
digital simulations provide detailed measurements on shorter time scales and smaller networks. 

We are pursuing a number of potential enhancements to our J J Neuron circuit. Along with existing tunable circuit 
parameters that affect the threshold, refractory period, frequency, size and shape of APs, more flexibility is possible 
by considering circuits with more junctions representing additional ion channels. To the extent that bursting and 
other exotic neuron behavior is caused by relatively slow ion channels, adding a large capacitively shunted junction 
(with smaller characteristic frequency) should allow JJ Neurons to display these behaviors. Plasticity of synapses 
based on previous history may also be possible. Connection strength can be dynamically adjusted with additional 
circuitry in the synapse. If the adjustments are dependent on the firing history of the neurons involved, a Hebbian 
learning process would be possible. 
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